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The in-medium behavior of the nucleon spectral density including self-energies is revisited within 
the framework of QCD sum rules. Special emphasis is given to the density dependence of four-quark 
condensates. A complete catalog of four-quark condensates is presented and relations among them 
are derived. Generic differences of such four-quark condensates occurring in QCD sum rules for light 
baryons and light vector mesons are discussed. 

1 Introduction 

A goal of contemporary hadron physics is to relate the confined quark and gluon degrees of freedom 
and parameters related to Quantum Chromodynamics (QCD) to the comprehensive hadronic spectrum. 
Lattice QCD and chiral effective field theory are suitable tools to accomplish this and other goals in 
exploring the structure of low-energy QCD and properties of hadrons. Another - though not so direct 
- but successful approach is given by QCD sum rules, originally formulated by Shifman, Vainshtein 
and Zakharov [HI] to describe masses of light vector mesons for example 10] ■ The method since then 
gained attention in numerous applications, e. g. to calculate masses and couplings of low-lying hadrons, 
magnetic moments, etc. (cf. e.g. flBS)- Its particular meaning is that numerous hadronic observables 
are directly linked to a set of fundamental QCD quantities, the condensates and moments of parton 
distributions. 

Hadrons are excitations from the ground state. Changes in this state are expected to reflect in a change 
of hadronic properties, especially in spectral functions and moments thereof related, e.g., to masses of 
hadrons. Measurements of "mass modifications" of hadrons in a finite temperature, strongly interacting 
medium or when situated inside nuclear matter, that means embedded in a bulk of protons and neu- 
trons and baryonic and mesonic resonances, then probe the QCD vacuum (cf. [6] for an experimental 
overview). The properties of the QCD ground state are mapped to and quantified by a number of con- 
densates, like intrinsic material constants, which partially carry information on symmetry features of the 
theory. Besides the chiral condensate (gg), up to mass dimension 6 the gluon condensate ^^G^), the 
mixed quark-gluon condensate {qgsCrGq), the triple gluon condensate (ggG^) and structures of the form 
{qTqqVq) contribute in vacuum. (F denotes all possible structures formed by Dirac and Gell-Mann ma- 
trices.) We emphasize here the specific role of the latter class of hitherto poorly known condensates, the 
four-quark condensates. 

Within the realm of hadron spectroscopy the explanation of the actual numerical value of the nucleon 
mass is crucial as ingredient for understanding macroscopic matter. The nucleon represents as carrier 
of mass the hard core of visible matter in the universe and thus is an important source for gravitation. 
Our investigations here are to be considered in line with previous investigations [0, S, for nucleons 
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inside cold nuclear matter, which are also discussed in and continuously explored in [|l2l. ll3L ll4l.ll5 



16l . llTL llSL Il9|] . The possible extension to further effects at finite temperatures lies beyond the present 



scope. An advantage of applying QCD sum rules is that effects of small finite baryon density n are 
described systematically by the change of condensates and the advent of new structures which are absent 
at zero density. For the nucleon an important dependence of self-energies on four-quark condensates 
was found. Comparisons with results of chiral effective field theory |120(], where nucleon self-energies 
show strong cancellation effects (i. e. they change with the same magnitude but have opposite signs) 
suggest that the relevant four-quark condensates should be weakly density dependent [7]. For the uj 
meson, however, we recently deduced in [21] evidence for a significant density dependence of a particular 
combination of four-quark condensates appearing there. This was based on the experimentally found shift 
of the spectral strength of uj to lower invariant masses inside nuclear matter, as measured by the CB-TAPS 
collaboration [|22l] . Therefore, we will spell out explicitly the four-quark condensates in nucleon sum rules 
in medium, which up to now are usually given in the factorization approximation or are determined by 
special models. So one can directly distinguish the four-quark condensate structures in sum rules for 
vector mesons and baryons. 

The work is organized as follows. In section 2 we review the operator product expansion for the 
nucleon and the hadronic model to write out the QCD sum rule equations. The four-quark condensates 
are discussed in section 3 where an exhaustive list is presented and relations among these condensates are 
given. Afterwards we present a numerical analysis of their influences in a nucleon sum rule evaluation 
and compare to other results (section 4). Conclusions can be found in section 5. In the appendices, an ex- 
planation of the calculation of an OPE and remarks on four-quark condensate relations are supplemented. 



2 QCD Sum Rules for the Nucleon 
2.1 Current- Current Correlator 

QCD sum rules (QSR) link hadronic observables and expectation values of quark and gluon operators. 
This allows to determine properties of the low-lying hadronic excitations of the QCD ground state |\1'). It 
relies on the concept of semi-local quark-hadron duality applied to the time-ordered correlation function 

which describes the propagation of a hadron created from the vacuum. On one side it can be calculated 
in terms of quarks and gluons via an operator product expansion (OPE) for large space-like momenta 
q^. This introduces Wilson coefficients multiplied by local normal ordered expectation values of quark 
and gluon fields - the QCD condensates. Thereby the hadron is assigned an interpolating field rj which 
resembles the right quantum numbers and is built from the fundamental degrees of freedom of QCD. On 
the other side, the interpolating field couples to the hadron excited from the ground state and the corre- 
lation function can be related solely to the hadronic properties for > 0. By means of analyticity of 
the correlator n(g), dispersion relations equate both approaches. This leads to the celebrated QCD sum 
rules ill]. Analysing transformed dispersion relations in a suitable range of momenta allows a determi- 
nation of hadronic properties. The generalization to nuclear matter with non- vanishing temperature or 
chemical potentials relies on Gibbs averaged expectation values |. . .| \1') instead of vacuum expecta- 
tion values (0 |. . .| 0) ||23l]. In what follows we focus on the nucleon, calculate the OPE and discuss the 
hadronic side; the sum rule finishes this section. 
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2.2 Interpolating Fields 

Following an argument of loffe one can write down two interpolating fields representing 

the nucleon with the corresponding quantum numbers I{J^) = ^(l^)^ ^'^''^[UaC'JiMUbl^^j'^dc and 
^abc^^T(j^^^^^j^^^fiu^^ » when restricting to fields that contain no derivatives and couple to spin i onlylil 
Extended forms of the nucleon current may include derivatives [0, or make use of tensor inter- 



polating fields [|28L 129] (also used to extrapolate the vacuum nucleon mass via QCD sum rules [30] to 



(unphysical) larger values obtained on the lattice, comparable to similar efforts within chiral perturba- 



tion theory Oil] ). The complications in nucleon sum rules can further be dealt with when taking into 
consideration the coupling of positive and negative parity states to the nucleon interpolating field [|32i] . 

In this work, our structures are always written for the proton; by exchanging u and d the neutron 
is obtained (even the neutron-proton mass difference has been analyzed in this frameworkyj]). As 
interpolating fields a Fierz rearranged and thus simplified linear combination is widely used ITI] 

ri{x) = 2e"^^ {t[Mj(x)C754(x)]M,(x) + [ulix)Cdbix)]-f5Ucix)} , (2) 

which in the above basis reads 

V{x) = ^e'^"' {(1 - t)[ul{x)C-fM^)HYdc{x) + (1 + t)[u];{x)Ca,Mx)]-f,a^''d,{x)} . (3) 

Both currents, r] and r/, are related by Fierz transformations whereby in such a straightforward calcula- 
tion the remaining difference vanishes for symmetry reasons (analog to the exclusion of Dirac structures 
in [|24] when constructing all possible nucleon fields). The consequence of these two equivalent repre- 
sentations ^ and (|3]) is that two different forms of the OPE arise. On the level of four-quark condensates 
the identity is not obvious and leads to relations between different four-quark structures. There appear 
constraints on pure flavor four-quark condensates which can be understood also without connection to an 
OPE (algebraic relations on the operator level). For mixed structures such relations do not follow. This 
will be discussed in section 3. 

Our subsequent equations will be given for the ansatz ^ with arbitrary mixing parameter t. In nu- 
cleon sum rule calculations the particular choice of the field with t = —1, the so-called loffe interpolating 
field, is preferred for reasons of applicability of the method and numerical stability of the evaluation pro- 
cedure (cf. also [34] for a discussion of an optimal nucleon interpolating field; another choice of t would 



emphasize the negative-parity state in the sum rule n35\\ ) 



2.3 Operator Product Expansion 

Using Wilson's OPE the correlation function ([T]) can be represented asymptotically as a series of Wilson 
coefficients multiplied by expectation values of quark and gluon operators, the condensates. As outlined 
in appendix 16. 1[ where also notations are summarized, these coefficients can be calculated considering 
quark propagation in a gluon background field which further simplifies in the Fock-Schwinger gauge. 
Application of Wick's theorem to ^ introduces the normal ordered expectation values, which projected 
on color singlets, Dirac and Lorentz scalars and restricted by the demand for time and parity reversal 
invariance in cold nuclear matter leads to an expansion into local condensates. 

The OPE for 11 (x) and the Fourier transform in ([T]) are important steps towards the sum rule formula- 
tion. Still the correlator can be decomposed into invariant functions. Lorentz invariance and the requested 
symmetry with respect to time/parity reversal allow the decomposition 

n(g) = n,(g2, qv) + U^(q^, qv)i + Il,{q\ qv)i , (4) 

'Concerning conventions on metric, Dirac and Gell-Mann matrices, charge conjugation matrix C etc. we follow ll25h . 
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where v is the four- velocity vector of the medium. The three invariant functions which accordingly yield 
three sum rule equations can be projected out by appropriate Dirac traces 

^.(g^g^;)=^Tr(^(g)) , (5) 
M<l^^v) = ^^^^^^ \ ^^^^^^ {^^^Tr (^n(g)) - MTr (^n(g))} , (6) 

=^-I-^^^l^^^ {g'Tr (^n(g)) - MTr(^n(g))} (7) 

and are furthermore decomposed into even (e) and odd (o) parts w.r.t. qv 

qv) = (qvf) + (qvf) . (8) 

For the nucleon interpolating field this leads to 

{qvf) = + H-q') {m) + ^ In(-g^) {qgs{aG)q) 

+ dk^^ ((^-(-^)^^/-^> + \ ^q9s{-G)q)^ , (9) 



(10) 



2 \ c^{qq) +^{qiq) \ „ , (H) 



- (l + ^) ((^-^(-^)^^/-^> + ^ {q9sl^i-G)q)^ , (12) 

+ ((gy(-/^)^g/^^^> + ^ (g^.y(-G)g)) , (13) 



^{qiq)l , (14) 



1 1 

where the Cj(t)'s, being polynomials of the mixing parameter t, are written out below in section [23] 
in the final sum rules. Numerical values for condensates are collected in section 14.21 where sum rules 
are numerically analyzed. The contributions from four-quark condensates are written here as the usual 
factorized result denoted by {. . .j^^'^; full expressions which replace and overcome this simplification 
are the focus of section |3] (see especially Eqs. (|501)-(|521) below). Note that, in contrast to the OPE for 
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n for light vector mesons with conserved currents = ^{ifj^u ± c^y^rf) (for u, p mesons), four-quark 
condensates enter already without a factor (the strong coupling) and the chiral condensate (qq) does 
not appear in a renormalization invariant combination with the quark mass. 



2.4 Dispersion Relations 

The representation of the correlation function through Wilson coefficients and condensates valid for large 
Euclidean momenta < is by the provided analiticity of 11 (g) related to the spectral density integrated 
over real values of the energy go- This is reflected in the fixed-g dispersion relation of the form (up to 
subtractions not displayed here and written now in nuclear matter rest frame only with go as argument) 

- / duj ^ = n go , (15) 

where g is held fixed and the spectral density p{uj) = ;^ An(c<;) enters as the discontinuity of the correlator 
n on the real axis 

AU{u) = ^ lim [U{uj + ie) - Ii{uj - ie)] . (16) 

Although dispersion relations could require polynomial subtractions enforcing convergence, such finite 
polynomials vanish under Borel transformation B : fiq"^) /(A^^) and need not be considered here. 
The correlation function is decomposed into even and odd parts using (HI) for f = (1, 0, 0, 0) defining the 
rest frame of nuclear matter, with 

n^(go^) ^ I (n(go) + n(-go)) = - du;^^^^ , (17) 

W{ql) ^ (n(go) - n(-go)) = - / du-^ , (18) 
2go uj^ - g^ 

which, given the integral representation (fTSl) . become functions of q^. The starting point of our sum rule 
analysis is the special combination 

- r^du;{uj-E)^^ = U%ql)-EU^{ql), (19) 

^ J-oo ^ — % 

where the l.h.s. encodes hadronic properties, while the r.h.s. is subject to the OPE (l9l)- (fT4l) . Such a com- 
bination has been proposed in f?] with the motivation to suppress anti-nucleon contributions effectively 
by a suitable choice of the quantity E. Having in mind the usual decomposition "resonance -i- contin- 
uum", we split the l.h.s. of the integral into an anti-nucleon continuum, — oo < uj < u;_, anti-nucleon, 
ci;_ < u; < 0, nucleon, < u; < and nucleon continuum, ujj^ < lo < oo, and choose 

E = —n and E = ^ -V^ j— ^ , (20) 
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where the latter similarly defines the first moment of the Borel weighted spectral density for the positive 
energy excitation]^ This delivers the Borel transformed sum rule 



(E-E)- dujAU{uj)e-'^'/^' =W{M^)-- dujum^^{u)e-'^'/-^' 



(21) 



The continuum contributions np'°j.(co') = AIl{uj) =i= An(— cj) are arranged on the r.h.s. with the reason- 
ing of employing the semi-local duality hypothesis: For the integrated continuum the same asymptotic 
behavior is assumed for the correlation functions of hadronic and quark/gluon degrees of freedom in the 
limit of large energies. These integrals are then extended to the respective continuum thresholds uj±. 
Typically only the logarithmic terms in 11 provide discontinuities which enter the continuum integrals. 
To summarize, Eq. (|2TI) exhibits the typical structure of QCD sum rules: the hadronic properties on the 
l.h.s., i.e. the low-lying nucleon spectral function, are thought to be given by the operator product rep- 
resentation of n including condensates on the r.h.s. The last term on the r.h.s. accounts for asymmetric 
continuum thresholds, i.e. co^ ^ —uj+, and can be estimated by semi-local quark-hadron duality. 

It should be emphasized that the given sum rule is for a certain, weighted moment of a part of the 
nucleon spectral function. Without further assumptions, local properties of An(c<j) cannot be deduced. 
Note also that in this form the anti-nucleon enters inevitably the sum rule. The reasoning behind the 
choice of (fT9l ) with (l20l) is that in mean field approximation, where self-energy contributions in the 
propagator are real and energy-momentum independent (cf. also [7J), the pole contribution of the nucleon 
propagator G{q) = {(^ — — H)^^ can be written as 

Gia)- ^ i + Mlj-i^. 

^^'^ - l-S,(,o-^.)(go-i^-) • ^''^ 
Pauli corrections to positive-energy baryons and propagation of holes in the Fermi sea give rise to an 
additional piece G'd(?) ~ ©(I^fI ~ 1^1) Il36ll vanishing for nucleon momenta q above the Fermi surface 
\(If\ considered here. The self-energy S is decomposed into invariant structures S = + '^qi + 
1371] (for mean field = 0), where one introduces scalar S., = — Mjy and vector self- 
energies Tjy, which are related to the decomposition above via M^^ = ^{^^^^ and = yrs" 
In the rest frame of nuclear matter the energy of the nucleon is E+, correspondingly E^ that of the 
antinucleon excitation, where E± = S„ ± ^/q"^ + . Since the sum rule explicitly depends on the 
nucleon momentum however the self-energy S as well as invariant structures E, and derived quantities 
acquire now a momentum dependence and become functions of the Lorentz invariants q^, qv and 
extending mean field theory towards the relativistic Hartree-Fock approximation [36]. Eq. (|22l) is giving 
rise to a discontinuity AG'(go) = ^ linie^o(G'(go + ^e) — ^(^0 — ^e)) with a simple pole structure 

AG(go) = ,^^f^E^'t~J^' (^(^0 - E_) - 5(go - E^)) , (23) 

where the general expression, Eq. (|20|) . identifies E with the anti-nucleon pole energy E^ for all 3 Dirac 
structures (analogously, E is identified with E+). Then the l.h.s. of the sum rule ((2T)) reads 



{E+ - EJ)- / dujM{{uj)e-^^'^^ = + ^^*N - ^S,)e-^+/^' . (24) 



TT 



1 — 



^In general, the Dirac structure of AH would require definitions Ei, Ei for the distinct invariant functions (i = s, q, v) 
of the decomposition (|4]l. In the case considered here we assume that these weighted moments coincide with Eg q y = 
E (analogously E). Also uj± are simplified to be common for s,q,v parts. In the shown Borel transformed equations, 
decomposed terms are symbolically rearranged to full Dirac structures. 
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Here, \n enters through the transition from Eq. ^ to the nuclon propagator via insertion of a com- 
plete basis, retaining only the nucleon state \q, s) with the relation to the nucleon bispinor Xiyu{p, s) = 

(\E'|?7(0)|g, s) and can be combined into an effective coupling = (XZs^- More general one can inter- 
pret Eq. (|24|) as parametrization of the l.h.s. of (|2T1) . where integrated information of AH is mapped onto 
the quantities M^, S^, and E±, which are subject of our further analysis, by virtue of Eqs. ©-dV]) 



\*^T.,e-''ll^" = {E 



1 

Jo 
Jo 

Jo 



(25) 
(26) 
(27) 



Due to the supposed pole structure in (|22)) the self-energy components are related to E± (or more general 
to E and E) and the relations from distinct Dirac structures are coupled equations. The given general 
spectral integrals however not yet relate the unknown quantities, so that our numerical results presented 
here are not completely independent of the given nucleon propagator ansatz. These relations highlight 
also the dependence on the Borel mass M. which one gets rid of by averaging in an appropriate Borel 
window. 

In [IbsIi . it has been pointed out, that 11 also contains chiral logarithms, e.g. m^^log rr^, which, 
however, do not appear in the chiral perturbation theory expression for M^. It was argued 11391. 14011 that 
low-lying continuum like txN excitations around cancel such unwanted pieces. In this respect, the 
parameters A/^, S^, S^, in (I25l)-(l27l) are hardly to be identified with pure nucleon pole characteristics, 
but should be considered as measure of integrated strength of nucleon like excitations in a given interval. 
Moreover, many hadronic models point to a quite distributed strength or even multi-peak structures (e.g. 
[|4lll ). The importance of an explicit inclusion of scattering contributions in the interval . . . oj+ has been 
demonstrated in [1421 |43L 14411 for finite temperature effects on the in-medium nucleon or in [45] for the 
in-medium when trying to isolate the pure pole contribution. In vacuum QCD sum rules for baryons, 
e. g. the nucleon, improvement of the continuum treatment is achieved by the inclusion of negative-parity 
states, which are equally described by a given correlation function as the corresponding positive-parity 
states [1321 |35L |46L |47L 14811 . Resorting to integrated strength distributions avoids these problems, but loses 
the tight relation to simple pole parameters. 



2.5 Sum Rule Equations 

Eq. (|2TI) is the sum rule we are going to evaluate with respect to the above motivated identifications. 
Inserting the decomposition (HI) with (l5l)-(fT4l) we arrive at the three coupled sum rule equations 

A^M^ + A^W + A3 , (28) 
B^M^ + B2M^ + 53 + B^/M^ , (29) 

CiM^ + + C3 , (30) 



'{El-q^)/M'^ 



1 



with coefficients 



A, 
Bo 
B2 



Cl 



— {gsqcrGq) 



2 UqiDQiDoq) + ^ {gsqcrGqU - {ci (qq) {qiq)}l 



37r2 
C4 



eff ' 



2567r4 
C4E0 



E- (q^q) 



+ 



247r2 
C4E0 I a. 



bc^Eo / I \ C4E0 /tts^2 



2567r2 \ TT 



11527r2 \ TT 

'^4^^ / t - n \ 
— ^ {qhDoq) 



288772 \ TT 



727r2 
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eflf 



^4 



C4 



El (gtg) 



127r2 

= ^,B,E^ W^Doq) - ^E_ (^[ivOr + (vG) 



2887r2 



TT 



C5E0 



c. 



C4 

27r' 



487r2 

■q' (^{qHDoiDoq) + ^ (^7.gVGg>^ ^l^A^ ^^"^^^T 



[gsq^f^Gq) 



eff 



and factors ii^, emerging from continuum contributions, witli tlie definition sq = — q'^. 



En 



^1 



E, 



So 



+ 



M2 2M4 



(31) 



-So/M2 



(32) 



and the asymmetric continuum threshold integral in Eq. (I2TI) neglected. The list (1311) is exhaustive for all 
condensates up to mass dimension 5 in the limit of vanishing quark masses. The coefficients q denote 
general structures due to the mixing of interpolating fields according to ^ obeying 



Cl = 7r - 2t - 5, 

C2 



-t^ + 1, 



C3 = 2t-t- 1, 

C4 = St^ + 2t + 5, 
C5 = 7t^ + lot + 7. 



(33) 
(34) 
(35) 
(36) 
(37) 



3 Four-Quark Condensates 

Formally, four-quark condensates are QCD ground state expectation values of hermitian products of four 
quark operators which are to be Dirac and Lorentz scalars, color singlets and are to be invariant under 
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time and parity reversal. Thereby we restrict ourselves to equilibrated cold nuclear matteitl but do not 
impose isospin symmetry from the very beginning in view of further applications, such as the proton- 
neutron mass difference in asymmetric cold nuclear matter (e.g. [17]). With the following discussion of 
independent four-quark condensates for arbitrary numbers of flavors we allow for the inclusion of strange 
quark contributions as well. Physically, the four-quark condensates quantify the correlated production of 
two quark-antiquark pairs in the physical vacuum. In contrast to the square of the two-quark condensate, 
which accounts for uncorrelated production of two of these pairs, the four-quark condensates are a mea- 
sure of the correlation and thus evidence the complexity of the QCD ground state. Especially, deviations 
from factorization, the approximation of unknown four-quark condensates in terms of the squared chiral 
condensate justified in the large Nc limit (cf. also [49]), represent effects of these more involved correla- 
tions. In this section the classification of four-quark condensates, in the light quark sector, is performed 
in some detail. 



3.1 Projection and Classification 

The projections onto Dirac, Lorentz and color structures lead to all possible in-medium four-quark con- 
densates just as for the example of the non-local two-quark expectation value in appendix 16. II However 
the situation is even simpler since we are only interested in the mass dimension 6 condensates, so deriva- 
tives are not required and all operators in four-quark expectation values are to be taken at x = 0. 

Using the Clifford bases G {1, 7^, (T/,<,,, ^757^5 75} and O™ G {1, 7^", a^'<'' , ^5!^, 75} which ful- 
fill Tr [OkO"^) = 45^ one can project out the Dirac indices of products of four arbitrary quark operators 

/ V 16 

gi"'g2'^g3''g4' = tt^ E Uio^q^.^O^qi) & Oi . (38) 

V e f g h / 16 V / f,e h,g 

Note, here and elsewhere, Dirac indices, if explicitly shown, are attached below the concerned objects. 
From (l38l) there are 25 combinatorial Lorentz structures which have to be projected on condensates to 
obey Lorentz invariance (using the four-velocity t;^), time/parity reversal and hermiticity. For each of the 
remaining 5 (10) Lorentz scalars in vacuum (medium) two possible color singlet combinations can be 
formed using contractions with the unity element and the generators of SU{Nc = 3). Thus one obtains 
the projection formula 

^i'g2^3'94 = I {<liq2qm) laa'lw + ^ (^i A^ga^s A^g4) Af,,Af,, . (39) 

Especially, in the calculation of an operator product expansion for baryons the color condensate structures 
naturally arise from the product eabc^a'b'd = ^abc^dVc = ^aa'^bb' — ^ab'K'b, hence there the four-quark 
condensates generally appear in linear combinations of color structures in the form 

eabcea'b'c g?'g2^3'?4 = ^ |(^i?2g3g4) - ^ (giA'^g2g3A'^g4) I • (40) 

This would imply two condensate structures for each Lorentz scalar term; however, for expectation 
values with just one flavor (pure flavor four-quark condensates) these structures are not independent. 
Combining Fierz rearrangement of the Dirac contractions of pure four-quark operators with the rear- 
rangement of the color structures, one derives the transformation equation 

(uOkX^uuO^X^u) = -I (uOkUuO^u) - ^Tr (OkOnO^O"") {uO^uWu) , (41) 

3 8 



■'The catalog can be extended to non-equilibrated systems lifting the demand for time reversal symmetry or to systems at 
finite temperature and vanishing chemical potential where charge conjugation provides a good symmetry. 
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which relates the two different color combinations. This transformation can be brought in matrix form 
y = Ax with 



/ {qX^qqX^q) \ 

{q-faX^qqj°'X^q) 
(q^X'^qqi/iX^q) jv^ 
{qa^pX^qa'^^X^'q) 
{qa^pX^qqa^'X\) g'^v^s/v^ 
{q'l5laX'^qq'y5l°'X^q) 
(qjniX'^qq-f^^X'^q) jv^ 

{ci'jX^qqX^q) 



X 



( (qqqq) \ 

{Qlaqqi'^q) 
{qi/iqqi/!q) jv^ 
{qaapqqa'^f^q) 

{qibiaqqibTq) 
{qi^iqqi^iq) h"^ 
{qimim) 
{qiqqq) 
\{qibTqq(^'^^q) iea^p^sv^ 12) 



-7/6 


-1/2 





-1/4 





1/2 





-1/2 





\ 


-2 


1/3 











1 





2 








-1/2 


1/2 


-5/3 


-1/4 


1 


1/2 


-1 


1/2 








-6 








1/3 











-6 








-3/2 


-1/2 


2 


1/4 


-2/3 


1/2 


-2 


-3/2 








2 


1 











1/3 





-2 








1/2 


1/2 


-1 


1/4 


-1 


1/2 


-5/3 


-1/2 








-1/2 


1/2 





-1/4 





-1/2 





-7/6 
































-5/3 


—i 


























3i 


i/s) 



(43) 



We emphasize that the inverse transformation A^^ exists. However, structures for baryon sum rules 
typically are combinations of two color contractions, dictated by Eq. (|40|) . which form components of the 
vector 

z=^(^x-y^ =Bx = ^(^A-'-^l^y, (44) 

The matrix B has the fivefold eigenvalues and 2, and the corresponding eigenspaces both have dimen- 
sion 5, especially the kernel of 13 spanned by the eigenvectors to eigenvalue 0. The fact that the kernel 
contains more than the null vector implies that B has no inverse. The transformation of this equation into 
the basis of eigenvectors yields a new vector z' where 5 elements are to be zero. Written in components 
of z these relations are 



Z2 + zg = 0, (45) 

4zi -2z2-z^ = 0, (46) 

2zi - Z4 + 2^8 = , (47) 

zi- z^~ Z5 + zj = 0, (48) 

zg - izio = . (49) 
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The first three conditions occur already in the vacuum set, the latter two ones are additional in the medium 
case. Of course, the conditions can be written differently, e.g., the second and third line may be conve- 
niently combined to zi — Z2 — Zs = for applications. An alternative derivation of these relations is 
presented in appendix l6.2[ 

The relations (l45l)- (|49l) have two important consequences: firstly, they allow to simplify pure flavor 
four-quark condensates in baryon sum rules; secondly, since Eq. (|44l) can not be inverted, they forbid a 
direct translation from pure flavor four-quark condensates in baryon sum rules at the order a° to those 
which occur e.g. in sum rules for light vector mesons in the order a]. 

3.2 Four-Quark Condensates in the Nucleon QCD Sum Rule 

We have now provided all prerequisites to specify the four-quark condensates which occur in the sum 
rule (I2TI) . The full expressions for the four-quark condensates in the order abbreviated in (|3T1) so far 
symbolically, are 

{ci {qq) mnls = l^abcea'b'c (-2c2 (m"Vw"w''w') + C6 {u^'iu''d''d') - 3C2 (tx^'m^J^V^') 

+C7 (u'''^^^^u''d''ax^d'e^^^%)) , (50) 



ci {qq? + ^ (^^?)'}'^ = ^abcea'b'c (2C9 lu'^' ^^u'^u'' ^^u') - 2c9 (m'^Vw^w'VwV^') 



- 9c2 {u^'u^d^'d^^ + ^C2 (u^'a^xu^d^'a^^d^^ - 9c2 (u^' -i^u^d!"' -id^ 
+ Clo (u'''-frU''d^'Yd^) - 2c9 l^'^iu'^^'id^jv'^ 
+C8 (n'^'757rM'^rf'''757"^^') - 4t ^xf' ^^u^^' ^r^d'' jv")^ , (51) 

- 2t (^u'*'757rM"^^'''757^M^) + 8t (u'''-f5iu''u'''-f5iu^/v^^ 

- C9 (u'''-fru''d^'Yd^) + 4c9 (u'''iu''d^'id^/v^^ 

-2t (^u^'^^^^u^d^'^^Yd') + 8t (u^'^^iu^d^'^^^d'/v^^^ . (52) 
Here, additional polynomials which express the mixing of interpolating fields are 

CQ = t^ -2t + 1, (53) 
C7 = - t, (54) 
cs = 9f + lot + 9, (55) 



C9 = r + 1, (56) 
cio = llt^ + 6t + 11. (57) 

These expressions extend the non-factorized four-quark condensates for the nucleon in vacuum listed 



in[|42L|5011 
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3.3 Factorization and Parametrization of Four-Quark Condensates 



Up to now we have introduced all possible four-quark condensates in the light quark sector and written 
out explicitly the structures which appear in the nucleon sum rule. In such a way the sum rule equations 
of the type employed, e.g. in [10], are equipped with complete four-quark condensates. We evaluate 
now the sum rule equations with the focus on these particular combinations of four-quark condensates. 
So we are faced with the common problem of the poor knowledge of four-quark condensates. Usually 
assuming the vacuum saturation hypothesis or resorting to the large Nc limit the four-quark condensates 
are factorized into products of condensates with two quark operators. The factorization of four-quark 
condensates allows to set the proper units, however its reliability is a matter of debate. For instance. 



Din state that the four-quark condensates in the nucleon sum rule are the expectation value of a chirally 
invariant operator, while (gg)^ is not invariant and thus a substitution by the factorized form would be 
inconsistent with the chiral perturbation theory expression for the nucleon self-energy. The four-quark 
condensates breaking chiral symmetry might have a meaningful connection to the chiral condensate but 



for the chirally invariant structures such a closer relation to {qq) is not clear [|52|] . 

Moreover, for nucleon sum rules at finite temperature T (and vanishing chemical potential) it was 
argued in [42] that the four-quark condensates are T independent quite different from the behavior of 
(gg)^ which is why a naive factorization would lead to artificial temperature effects in the nucleon mass. 

For numerical purposes it is convenient to correct the values deduced from factorization by factors k 
and examine the effect of these correction factors on predictions from QCD sum rules. In this section the 
four-quark condensates classified so far in general are spelled out and the parametrization with a set of 
quantities k is defined. In doing so one includes a density dependent factor K{n) in the factorized result 

(g/iriCig/ig/2r2C2g/2) = {qfi^i'Ciqfiqf2^2'C2qf2)i^^ , (58) 

where k and the following parametrization depend on the specific condensate structure. In linear density 
approximation this product ansatz obtains contributions both from the expansion K{n) = k'^'^'> + K^^^n with 
K^^^ = and from the linearized, factorized four-quark condensate expression (g/iriCig/ig/2r2C2g/2)fac 
a + hn. If = l,thenKW = recovers the usual factorization, which means the four-quark condensate 
behaves like the product of two two-quark condensates; k^^'> > represents a stronger density depen- 
dence with respect to the factorization and vice versa. Inserting both expansions one can also describe 
the total density dependence of the condensates by the combination n^'^^ = n^'^^ + f k^^), 

(g/iriCig/ig/2r2C2g/2) = aK(°) + bn^^'^n (59) 



such that for k^'^'^ = the condensate is (in first order) independent of density. For condensates with 
vanishing a or 6 in factorization we choose a = {qq)1^^ and h = {qq)^^^^ a^/niq as scale to study deviations 
from zero and denote these instances by R. The classification of possible four-quark condensates is 
collected together with the specific k parametrization in Tabs.[i|and[2l 



Non-Flavor Mixing Case 

The condensates which contain only one flavor are listed in Tab. [B From the demand for parity and 
time reversal invariance only 5(10) Dirac and Lorentz scalar four quark operators remained in vacuum 
(medium). Further, these structures carry color indices and must be projected on colorless objects for 
which there are two ways. However, since the same flavors occur, both color combinations can be alter- 
natively rearranged via Fierz transformation. Hence, there are only 5 (10) independent k parameter sets 
in the Tab. [T} although both color alternatives are listed. The parameter sets with indices 1, 2 are related 
by the transformation (|43l) . 
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Flavor Mixing Case 

Here the condensates containing two quark operator pairs are distinguished by flavor. The numbering is 
as for the pure flavor structures. However, the conversion of the two color contractions is not possible 
due to different flavors. Compared to the non-flavor mixing case the missing exchange symmetry of qq 
contractions due to different flavors allows additional placements of Dirac matrices and thus leads to 
4 additional condensate structures in medium (see Tab. [2l). Therefore, 10 (24) flavor-mixed four quark 
condensates and thus k parameter pairs appear in vacuum (medium). 

Hence, there exist in medium {vacuum} for n j flavors without flavor symmetry taken into account 
2nf{6nf — 1) {5nj} independent four-quark condensates being Lorentz invariant expectation values of 
hermitian products of four quark operators constrained by time and parity reversal invariance. Symmetry 
under flavor rotation reduces these numbers to 20 {10}, respectively. Finally note that these are also the 
numbers of necessary k^'^^ parameters. Since the four-quark condensates in operator product expansions 
obtained from the medium projections in the limit of vanishing baryon density n should coincide with the 
vacuum result, this leads by contraction of vacuum and medium projections of four-quark condensates 
to the relations /t™^, ^, = which have already been included in Tabs. [Hand [21 Further, Lorentz 

projections which exist only in medium imply no new k™'^ parameters and so the number of k^^'^ in 
medium reduces consistently to the number of k™'^ and four-quark condensates in vacuum. 
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Indices Full condensate 



Parametrized Factorization 
in Linear Density Approximation 



Is 
Iv 
IV 
It 
It' 
la 
la' 

Ip 
Ivs 
lat 



(uuuu) 

{ui^uuu) 



12 

- (^K 
2 U'^ 



vac 

s 



vac 
Iv 



- K, 



vac 
2 

vac 
2 

vac 
2 



It 

vac 
It 



- (- 

12 V 4 

-a 



vac 
a 



i^vac 
4'^la 



K 



Ip 



im) 
{qq) 

(W)vac 

(W)vac 
(^?)vac 



"Is '"0 



"It 



2s 

2v 

2v' 

2t 

2t' 

2a 

2a' 

2p 

2vs 

2at 



(m7„A^ 
(ui^X 



^uwy^'X^u) 



(uGapX^UUG^'I^X^u) 

(uGapX'^uua'i'^ X'^u) g^v^vs/' 



{u'^blaX 
{u'-f^iX- 

('U75A^ 



(u'^X^uuX^u) 



{uj5j^X^UUGXnX^u) e'^^^'^v^ 



vac 
2s 



vac 
2v 



4/1 vac 



vac 
2 

vac 
2 



9 V4"'2v 
16 / i^vac 

"I (i'^ 

16 f i^vac 

9 U'^ 

4 /'^vac 
9 V'^2p 



^2t 



vac 
2t 



^2a 

vac 
2a 



(^?)vac 
(^?)vac 

(^?)vac 



+ 

+ /t^v''^) 



+ K. 



med 
2t 







mcd 
2a 







Z-raed 
'^2vs 



2at 



Table 1: Two complete sets (indices 1 and 2) of independent non-flavor-mixing four- 
quark condensates differing in color structure and their parametrization with k in strict 

linear density approximation = (w)vac '^N/?n^q)- The sets are related by a Fierz 
transformation. A similar table for flavor d instead u appears for an exhaustive list of 
four-quark condensates for the two-flavor case Uf = 2. 
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Indices Full condensate 



Parametrized Factorization 
in Linear Density Approximation 



3s {uudd) KIT {qq)L + '^sT^^e 

3v (uj^udrd) KIT + ^"v^< 

3V' {Uiudid)/V^ l^vac ^^^^2^^ + ^med^^ 

3t {ua^puda'^Pd) RY^^ {qq)l^ + kZ^'n^ 

3t' (ua^puda^'d) g^vs/v^ i^^f (^g)vac + 

3a' (t.75^^.rf75yt^> /v^ (^9)L + '^?a'^< 

3p (uj^udj^d) klT {qq)l^ + k^p'^n^ 

3vs {uiudd) (gg)^^^3n/2 

3at {u^^^McJx.d) e^^^v^ (gg)^^ 3n/2 

4s (uX^udX^d) kZ (gg)L + f^l'M 

4v' (uix^ud^x^d) /v^ Ikz {qq)Lc + '^■r''^^? 

4t {ua^pX^uda'^PX^d) k^ {qq)t, + k^M 

4t' {ua^fsX^uda-^'X^d) g^v^vg/v^ ^k^ + ^4t''< 

4a (iZ757aA^iiJ757°A^d> {qq)l, + k^^n^ 

4a' {u-f^X^'ud-f^iXH) jv-" \kz (gg)L + /^^a^^ 

4p iu^^x^uh^xH) «ir(?g)L + ^?pV 

4vs (iZ^A^i^JA^d) «?vs'*(g?) vac 3^/2 

4at (tZ757«A^iiJ(7A.A^d> e'^^-^^;^ (gg),^ 3n/2 

5vs (dVdiZii) Ac^^f (gg),^ 3n/2 

Sat ( J757Krffi<7A.^i> e"^"^^e (??)va. 3n/2 

6vs (JyA^dfZA^ti) k^^ {qq)^ 3n/2 

6at (d757KA^dfZa;,,A^ii> e'^^-^^;^ kZ^ (gg),^ 3n/2 



Table 2: A complete set of independent flavor-mixing four-quark condensates and their 
parametrization by k parameters in strict linear density approximation. Additional pa- 
rameters (indices 5 and 6) are required for structures which cannot be exchanged. 
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Insertion of these parametrization into the relevant sums of four-quark condensates (|50l) -(l52l) yields 
effective k parameters as linear combinations of the previously defined condensate- specific parameters. 
The sum rule is only sensitive to these effective combinations and can thus only reveal information on the 
behavior of specific linear combinations of four-quark condensates. Therefore in the sum rule analysis 
the three parameters k^'^'^, k'^^'^, k™^'^ describing the density dependence enter as 

{ci m m) le'fT = («:r' {m).^ \^ , (60) 

{ci {qq? + ^ {qiq?]\ = Ci (^T i^^)L + i^^)..c ^n) , (61) 
{q^qfY = (kT' iml.. —n] (62) 



J cS V "^q 



and are functions of the mixing angle t as well. However, we restrict this discussion to the limit of the 
loffe interpolating field t = —1. Note again, the k^^^ values are effective combinations representing 
the density dependence of the respective condensate lists (|50l)-(l52l) and thus negative k^"^^, a four-quark 
condensate behavior contrary to the factorization assumption, comprise cancellation effects within these 
condensate combinations. 

Density dependence of four- quark condensates from models 

It is instructive to derive values for the effective density dependence parameters k™'^'^. Expectation val- 
ues of four-quark operators in the nucleon were previously calculated in a perturbative chiral quark 



model [53] and taken into account in sum rule evaluations for the in-medium nucleon rtl6ll . (Correc 



tions to the factorization of four-quark condensates in nucleon sum rules have also been considered in the 



framework of the Nambu-Jona-Lasinio model in [|54l] . Lattice evaluations of four-quark operators in the 



nucleon are yet restricted to combinations which avoid the mixing with lower dimensional operators on 
the lattice [55], and provide not yet enough information to constrain Eqs. (l50l)- (|52l) .) The results in [53] 
can be translated to our k parameters. However, only such color combinations being significant in baryon 



sum rules are considered, see left column in Tab. [3l We note that the values given in 115 311 have to be cor- 
rected slightly in order to reach full consistency with the Fierz relations (|45l)-(|49l), which are an operator 
identity and thus must be fulfilled also for expectation values in the nucleon. An optimized minimally 
corrected set is found by the following procedure: minimize the relative deviation of all separate values 
compared to values delivered in the parametrization of [53] (this is in the order of 10 %, however with 
different possible adjustments); from these configurations choose the set with smallest sum of separate 
deviations (this deviation sum estimates to 40 % and different configurations are close to this value). 
The results from which the relevant density dependence for our condensate classification is obtained are 
collected in Tab.[3l our slight modifications of values in the original parametrization [53] are documented 
in Tabs . |4] and \5\ in appendix 16. 3[ 

The connection to our k parameters is derived as follows: Generally, in linear density approximation 
condensates behave like (\l/o \ A\ ^q) = (0 1^41 0)+n {N \ A\ N). If one compares our parametrized density 
dependent part of each four-quark condensate with the evaluation of nucleon matrix elements of four- 
quark operators in the combinations in Tab.[3]one obtains values for linear combinations of k parameters. 
The linear combinations refer to the two distinct color alternatives representing, as mentioned above, the 
typical color combination in baryon sum rules. These values can thus be applied to give also the required 
effective parameter^ (apart from the term ~ in"'' 'y^'y^u'^df'' axT.d!'') ^ e'^^'^^v^ not considered in [53] which 



Note some difference to the OPE part stated in eqations (87)-(89) of 01611 for the whole combination of the density depen 
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Mpan Niiclpon IVTatrix Plpmpnt 
(to be color contracted with Cabc^a'b'c') 


PCOM model \(ad) 1 

rx^V^^A iiiuuci [\yy/vacJ 








1.977 




0.432 




12.024 


(jQ^pUi U (J U 1 U U j U 






— 1.980 




—0 ^1 9 




2 Olfi 










in"-' u"-d!'' <f\ 


3.19 


a 7 a 




l^l"-' 4lfi(f>' 4(f>\ /7.2 


—0 73 




3.36 




1.11 




1.66 


('u'''75'^M''(?'75yc?'')^/t;2 


0.37 


(m"'75M'^(?'75C?'>^ 


-0.185 




-0.245 


(M'^'757KM"(i'VA.rf'')^e'^^<t;5 





Table 3: The combinations arranged as in the vector z oi four-quark expectation values 
obtained from the (partially modified) set taken from a perturbative chiral quark model 
calculation (PCQM) in [53] from which the characteristic density dependence of four- 
quark condensates, the value of is derived. Isospin symmetry N = i(p + n) of 
the nuclear matter ground state is assumed. The values in the pure flavor sector (upper 
part) are tuned to obey Fierz relations (l45])-( |49t on the accuracy level < 0.01 {qq)^^^- For 
three combinations no results are provided in i53n as indicated by "— ". 
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we had to neglect in the determination of ) 



K^"'^ = -0.25 , 



(63) 



Note that individual parameters are not small compared to these effective numbers indicating signif- 
icant cancellation effects in the density dependent parts of combined four-quark condensates. Moreover, 
for pure flavor four-quark condensates the ambiguity due to Fierz relations between operators does not 
allow to prefer a specific condensate type as dominating the four-quark condensates in the sum rule. For 
further attempts to gain estimates of four-quark condensates we refer the interested reader to 115 611 . 



Remark on Chiral Symmetry 

The chiral condensate is often considered an order parameter for the SU{nf )A chiral symmetry of QCD. 
Its change however might partially originate from virtual low-momentum pions and thus could not clearly 
signal partial restoration of chiral symmetry in matter 14011. The interpretation of four-quark condensates 
as order parameters for spontaneous break-down of chiral symmetry is an open issue. In ll52ll a specific 
combination of four-quark condensates arising from the difference between vector and axial-vector cor- 
relators is proposed as such an alternative parameter. This combination (shown to agree with vacuum 
factorization in the analysis of tau-lepton decay data ||57|]) is distinct from the above four-quark conden- 
sate lists in the nucleon channel as well as from the combination in the cu sum rule [21]. For instance, in 
vacuum nucleon QCD sum rules the four-quark condensate combination (the vacuum limit of Eq. (ISTI) 
with isospin symmetry being applied; ^ is the flavor vector) enters as sum of a chirally invariant part 



[2{2t^ + t + 2) (ipjf.iptpjf'tp) + {3f + 4t + 3) {'4'-f5l^^^i^l5Y^)] - ^ [color structures with A^] 



and a part which breaks this symmetry (pointed out in factorized form already in BSll ) 



— 7 [color structures with A^l 

4 ^ ^ 



In the preferred case t = —1 only the chirally invariant part survives and thus this remainder cannot 
be used as an order parameter. Additional insight into the change of four-quark condensates and their 
role as order parameters of spontaneous chiral symmetry breaking could be acquired from other hadronic 
channels, as the generalization of further baryon sum rules in vacuum Il58l.l59l.l60ll to the medium case, 
e.g. forthe A [|6lll62|l. 



4 Analysis 

4.1 Approximations 

In text book examples of QCD sum rules for light vector mesons (e.g. [jsi 0]) one usually considers mass 
equations and optimizes them for maximum flatness w.r.t. the Borel mass. This, however, includes often 
derivative sum rules and seems not to be appropriate in the case of fermions where the condensates are 
distributed over coupled sum rule equations for several invariant functions. Despite of this, equations 
for the self-energies can be formed dividing Eqs. (l28l) and (l30l) by (|29l) thus arriving at a generalization 

dent four-quark condensate contribution. Our equivalent OPE calculation utilizing the same nucleon four-quark expectation 
values (encoded in n^^f, k™'"^ as above) yields = (0.49^ 1 + 0.52^ + 0.57^^) -^ti, with p = Mnv. 
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of loffe's formula tlM for the nucleon vacuum mass. Approximated forms incorporating only lowest 
dimension condensates are sometimes used as estimates for in-medium nucleon self-energies |63i 64]: 

S„ = 647r2 (gtg) /{3M^) = 0.36 GeVn/riQ and = -M^ - {qq) /M^ = -0.37 GeVn/riQ at 
= 1 GeV\ 

Although to be confirmed by dedicated sum rule analysis, it is instructive to understand the impact of 
four-quark condensates at finite density from naive decoupled self-energy equations linearized in density. 
For fixed Borel mass ^A'^ = 1 GeV^, threshold sq = 2.5 GeV^ and condensates listed below, the self- 
energies become independent when a constant = —M^ is assumed; with k'^^ adjusted to yield the 
vacuum nucleon mass the self-energies are estimated as 

= (0.16 + 1.22^'^^) GeV— , (64) 

no 



-(0.32 + 0.11^^'='^ - 0.31k^"^<^) GeV— . (65) 



n 
no 

Indeed at small values of kp the impact of k^'^'^, k^'^'^ and k^'^'^ is as follows: The vector self-energy 
only depends on k^^'^, the scalar self-energy is effected by k™'^'^ and k'^"^, whereby a negative 
^med v^,orks equivalent to a positive value for k^^'^ and vice versa. Comparable effects in point out 
that a characteristic value of k™'^'^ is three times the corresponding absolute value of Whereas 
this qualitative estimate from Eqs. (|64|) and (|65l) is in line with the numerical analysis below for small 
densities n < O.Twq corresponding to Fermi momenta kp = (37r^r2/2)^/^ < 1.2 fm^^, the limit of 
constant four-quark condensates deviates from the widely excepted picture of cancelling vector and scalar 
self-energies which can be traced back to competing effects of higher order condensates. Since even in 
the small density limit for constant four-quark condensates the estimated ratio S^/S^ ~ 1 cannot be 
confirmed numerically, these estimates cannot substitute a numerical sum rule evaluation. 



4.2 Numerical Analysis 

In order to investigate numerically the importance of the three combinations of four-quark condensates 
entering the sum rule equations (l28l)- (|30l) at finite baryon density we perform an evaluation for fixed 
continuum threshold parameter sq = 2.5 GeV^ in a fixed Borel window Ai"^ = 0.8 . . . 1.4 GeV^. Since 
we are especially interested in medium modifications we use all sum rule equations although chiral- 
odd sum rule equations have been identified more reliable in the vacuum case [|65|] (however note that 



instanton contributions might change the relevance of particular sum rule equations [|66l 16711 ). From 
Eqs. (|28l) - (|30l) after transformation one unique left-hand side and the corresponding three r.h. sides are 
compared and their differences are minimized by a search for the optimum parameters S„, and 
with a logarithmic deviation measure [@, [ssl]. The condensates are estimated from various relations, 
e.g., the chiral condensate (qq) depends via partial conservation of the axial current (PCAC) on the pion 
decay constant and pion mass, the gluon condensate {{as/n)G'^) is determined from the QCD trace 
anomaly and further condensates can be expressed through moments of parton distribution functions. 
We use the values {qq) = -(0.245 GeV)^ + naN/i^rrig) with ctn = 45 MeV and nig = 5.5 MeV, 
(gsqaGq) = {qq) + 3.0 GeV^n with = 0.8 GeV^ {qiD^iDoq) + {gsqaGq) /8 = 0.3 GeV^n, 
{{as/7r)G^) = -2{{aJn){E^ - B^)) = -2[-0.5(0.33 GeV)^ + 0.325 GeVn)], {{aj7r)[{vGY + 
{vGf]) = -((a,/7r)(E2 + 52)) = 0.1 GeVn, (qUDoiDoq) + {gsq^ a G q) / 12 = (0.176 GeV)2n, 
{q^q) = 1.5n, (^qHDoq) = 0.18 GeVn and (gsq^crGq) = —0.33 GeV^n as employed and discussed 
in IS']. The values of possible k^^"^ parameters are given in Eqs. (|63l) . Values for are adjusted to 
reproduce the vacuum nucleon pole mass. 
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The results of these numerical evaluations for a nucleon on the Fermi surface Igp] = kp are sum- 
marized in Figs. [ill6l Fig. [T] shows the scalar and vector self-energies of the nucleon as a function of 
the Fermi momentum. The situation with four-quark condensate combinations (|60l) -(l62l) kept constant at 
their vacuum value (i. e. k^^^^ = k™'^*^ = k^'^^ = 0) is compared to the QCD sum rule evaluation with 
K parameters from Eqs. (l63l) . The results have the same qualitative behavior as self-energies determined 
from chiral effective field theory with realistic NN potentials poL 1681 ]. Figs. [2l [3] and |4] exhibit the im- 
pact of the 3 different four-quark condensate combinations: The vector self-energy is, in agreement with 
Eq. (l64l) . mainly determined by k™'^'^ especially for small densities (for positive values of even the 
qualitative form of the vector self-energy changes), k^"^^ has only small impact, and k,^^'^ does not effect 
Sv at all. The scalar self-energy, in contrast, is influenced by all 3 combinations, whereby the change of 
^mcd |g Qj^iy visible for Fermi momenta kp > 0.8 fm~^ as also suggested by Eq. (|65l) . Figs. [2] and [3] also 
reveal the opposed impact of n^^'^ versus k™*^*^. A variation of sq is not crucial (see Fig. [5]). The inclusion 
of anomalous dimension factors in the sum rule equations as in leads to a reduction of in the 

order of 20% but causes only minor changes in S^. Thereby the naive choice of the anomalous dimension 
from the factorized form of the four-quark condensates leaves space for improvement since it is known 
that four-quark condensates mix under renormalization [|69D. Our analysis concentrates on the impact 
of four-quark condensates, but also the variation of the density dependence of further condensates can 
change the result. For example, a large change of the density behavior of the genuine chiral condensate, 
as determined by the term, by factor 2 {0.5} leads to 8 % decrease {4 % increase} in the effective 
mass parameter at kp ~ 0.8 fm~^, while S„ is less sensitive. Correspondingly, the effective coupling 

is reduced by 10 % {enhanced by 5 %}. 

An improved weakly attractive cancellation pattern between (attraction) and S^, (repulsion), and 
thus agreement with chiral effective field theory [|68|] . can be achieved for a parameter set k^'^'^ = 1.2, 
K^'^'^ = —0.4, K^"^^ = 0.1 (see Fig. (6]). However such a fit would allow larger values of k^^'^ compensated 
by a larger magnitude of the negative value of k™'^'^ and vice versa. Note that in both ways the factoriza- 
tion limit K^^^ = 1 is violated by one or the other four-quark condensate combination. Such optimized 
K parameters, adjusted to results of [68], deviate noticeably from those in Eqs. (|63] ) deduced from ||53|] . 

Quantities characterizing the energy of an excitation with nucleon quantum numbers are and 
introduced in section 2.4. Since is negative, drops continously with increasing density achieving a 
value of about 540 MeV at nuclear saturation density (corresponding to kp ~ 1.35 fm^^) if extrapolated 
from the optimized fit in Fig. [6] The energy barely changes as function of kp. 

Considering the behavior of the effective coupling parameter in the cases of Figs. [2]l4] the maximum 
impact of k^^'^ {k^^'^} on A;^^ is 6 % {3 %} at kp ~ 0.8 fm'^. In the extreme case, k^""^ = 1 leads 
to a 40 % increase of A^. The variation of this coupling as a function of kp is in the order of 10 % 
in the optimized scenario. Generally, specific assumptions on the four-quark condensates can cause a 
decrease or an increase as well. This alternation of A^ has already been pointed out in [0], whereby their 
assumptions yield even a ±20 % change at nuclear density compared to the vacuum limit (cf. also n2m ). 
The vacuum limit of the calculated A^ agrees with the existing range of values (see [ 34|] for a compilation 
of results for the coupling strength of the nucleon excitation to the interpolating field in vacuum). 
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Figure 1: Nucleon vector and scalar self-energies as functions of the nucleon Fermi 
momentum kp = (37r^n/2)^/^. The sum rule result for constant four-quark condensates 
(QSR with constant fqc: k^^"^ = k^^'^ = k^^"^ = 0, solid curve) is compared to an 
evaluation with density dependent four-quark condensates as given in Eqs. (1631 ) (QSR 
with fqc from PCQM, dotted curves). The latter choice causes only minor differences 
in S„ and S^, for the scalar self-energy also because of competing impact of k^^^ and 
^med jjjg self-energies from chiral effective field theory i68ll (ChEFT, dashed curves) 
are shown as well but should be used as comparison only at small densities. 
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fcp [fm-1] 



Figure 2: The variation of nucleon self-energies for different assumptions of the density 
dependence of the four-quark condensates in Eq. (1281 ) parametrized by k™*''^ ; other four- 
quark condensate combinations are held constant. 
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Figure 4: The same as Fig.|2]but for a variation of k^'^^ {k^^^ = k^^*^ = 0). 
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Figure 5: The impact of different threshold parameters sq on the nucleon self-energies 
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Figure 6: QCD sum rule evaluations of nucleon self-energies with the parameter set 
K™'^'^ = 1.2, K™*^*^ = —0.4, K™'^'^ = 0.1 (dash-dotted curves) compare to chiral effective 
field theory P68l] with realistic NN forces as input. 



5 Conclusions 

Four-quark condensates have a surprisingly strong impact on conventional spectral QCD sum rules of 
light vector mesons. Unfortunately, four-quark condensates and their density dependencies are poorly 
known. One possibility is to consider a large set of hadronic observables and to try to constrain these 
para meters characterizing the QCD vacuum. Steps along this line of reasoning have been done, e. g. , in 



16211 . In order to accomplish a systematic approach, we present here a complete catalog of independent 
four-quark condensates for equilibrated symmetric or asymmetric nuclear matter. While the number of 
such condensates is fairly large already in the light quark sector, we point out that only special combina- 
tions enter the QCD sum rules. For the conventional nucleon QCD sum rule, three different combinations 
of four-quark condensates are identified. We note that the knowledge of these combinations (even the 
individual condensates entering) is not sufficient to convert them into the combination being specific for 
the spectral QCD sum rule for light vector mesons. In analyzing the set of independent four-quark con- 
densates we find also identities which must be fulfilled in a consistent treatment. Model calculations of 
four-quark condensates seem not to fulfill automatically these constraints. 

On the level of an exploratory study we show the impact of the three combinations of four-quark 
condensates on the vector and scalar self-energies of the nucleon. In cold nuclear matter at sufficiently 
low densities the density dependence of only one effective four-quark condensate combination is found 
to be important for the vector self-energy and the other two combinations dominate the scalar self-energy. 
Keeping in mind that the nucleon self-energy pieces per se are not proved to represent observables, one 
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is tempted to try an adjustment of these parameters to advanced nuclear matter calculations. While the 
overall pattern agrees fairly well (i. e. large and opposite scalar and vector self-energies) we can reproduce 
also the fine details on a quantitative level at low densities. Keeping the four-quark condensates frozen in 
to vacuum values or giving them a density dependence as suggested by a perturbative chiral quark model 
induce some quantitative modifications which may be considered as estimator of systematic uncertainties 
related to the four-quark sector. Furthermore, the special use of sum rules and interpolating current and 
details of the numerical evaluation procedure may prevent QCD sum rules for the nucleon as precision 
tool. The knowledge of this situation may be of relevance for approaches to the nuclear many-body 
problem which utilizes chiral dynamics and condensate-related features of the mean field. 

Finally, we remind that our study is restricted to cold nuclear matter. The extension towards finite 
temperature deserves separate investigations. 
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6 Appendix 

6.1 Operator Product Expansion 

For completeness and to clarify some technical details we recollect important steps of an OPE calculation. 
A convenient way to obtain this series is to calculate the Wilson coefficients in an external weak gluon 
field [70]. In the background field formalism the correlation function ([T]) is expanded according to Wick's 
theorem 11 (x) = nper(x) + n2q(a;) -|- n4q(a;) -|- . . . , where the full contractions are collected in the 
perturbative term Ilper and further terms n2q,4q,... denote the number of non-contracted quark operators. 
The latter terms give rise to non-local condensates containing the indicated number of quark operators. 



The use of Wick's theorem naturally introduces the normal ordering of operators ( \E' 



Ai---A 



n 



Ai---Anj, which will be assumed in all expectation values formed out of products of field operators. 

Under the presence of the gluon background field the quark propagator S'^ which appears in the terms 
in n(x) is modified, following from the solution of the Dirac equation in an external field in the Fock- 
Schwinger gauge for the gluon field. The corrections to the free quark operator appear in an expansion 
in the coupling Qs = v^vraj 

S'Jx) = \T [qa{x)m]\ ^) = 7^A^<^>^ + ^^5;^.(0)T„1^7^75 + • • • , (66) 

with the dual gluon field strength tensor G^^ = ^c^j^kaG'''^^ and color matrices T^, valid for massless 
quarks and inclusion of pure gluon condensates up to mass dimension 4. 

The Fock-Schwinger gauge is determined by {x — xo)f^Af^{x) = 0, and usually one chooses xq = 0. It 
allows to express partial derivatives of fields easily by covariant derivatives which matters when expand- 
ing non-local products of such operators. In general, results are gauge invariant, however technically 
fixing this gauge has enormous advantages in calculations of Wilson coefficients. Let us remark, that 
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although the term initially contains two uncontracted quark field operators, the expansion of the non- 
local expectation value into local condensates together with weak gluon fields resulting from modified 
quark propagators and the use of the equations of motion would induce further four-quark condensates at 
the order a^. 

The use of the quark propagator (|66l) leads to gluon insertions in the expectation values in 11 and thus 
to condensates of higher mass dimension. To obtain the condensates the expectation values are projected 
onto all possible Dirac, Lorentz and color scalars obeying symmetry w.r.t. time and parity reversal. This 
introduces all possible condensates up to the considered dimension, and having inserted the projections 
for the specific correlation function offers also the corresponding Wilson coefficients and therefore the 
OPE [8]. 

For example, the non-local diquark expectation value can be projected on color and Dirac structures 

{qaa{x)q,(,m = -^E^r (g(O)rg(x)) T^^ , (67) 

r 

where elements of the Clifford algebra F G {1, 7^, a^,y, ^757^*, 75}, are contracted over Lorentz indices, 
er = I for F = a^,y and er = 1 otherwise. A Taylor expansion of the quark operator at a; = in the 
Fock-Schwinger gauge 



q{x) = g(0) + xf^D^qiO) + -x^x'D^D,q{Q) + . . . (68) 



1 

leads to additional Lorentz structures, such that the local expansion of the non-local diquark term (|67l) up 
to mass dimension 5 in the expectation values taken at x = yields 



{qaa{x)qbp{Q)) = 



^vTap (^{qTq) + x^ (qTD^q) + ^x^x' {qTD^D.q)^ . (69) 



However, matrix elements (^(O)Fg(x)) with F G {cr^„, ^757^, 75} do not contribute due to the demand 
of time and parity reversal invariance and the multiplication with the symmetric Taylor expansion in 
X. Condensates with field derivatives can be transformed whereby a couple of manipulations using the 
equations of motion 

(z^-m)g = 0, g(z'^ + m) = 0, DfG^^ = g,Y.^YT\. (70) 

/ 

and the representation of the gluon tensor G^j,^ = TaG^^ 

G^, = -[D^,D,], and thus ^g^aG + I/)!^ = D^, = ^ (7^^ + ^7^) , (71) 

are exploited to yield simplifications in condensate projections. Terms which contain factors of the small 
quark mass are neglected in these considerations. 

Similar projections can be performed for structures which include gluonic parts from the propaga- 
tor (l66l) and lead to gluon condensates in npcr(a;) and are also carried out to find the linear combinations 
of four-quark condensates in 114^. Following this sketched line of manipulations, one arrives at Eqs. ([8])- 
dl). 
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6.2 Alternative Derivation of Pure-Flavor Four-Quark Condensate Interrelations 

The constraints between two different color structures of pure-flavor four-quark condensates have been 
presented in section [STI by analyzing the specific color structure transformation. For the typical baryon 
color combination of four-quark condensates the conversion matrix B (|44)) was derived with the decisive 
property that it cannot be inverted. In algebraic terms, the underlying system of linear equations is 
linearly dependent. This gave rise to the Fierz relations (l45l)-(|49l). If one is only interested in these 
relations, another direct way of derivation exists. Thereby one considers the "zero identity" 

^abc^a'b'c' q^'q^qb'qb (p^) ) = if {TCf = -{TC) or [CTf = -(Cf ), (72) 

e f g h c,g f,h 

which can be seen by a rearrangement of the product and renaming of indices (this is the analog discussion 
as for the choice of possible interpolating fields for the nucleon). Fierz transformation of this relations 
yields the basic formula 

^abc^a'b'c' (f'o^q-qVQnqh ^ {vOJ^CO'^^C^ = , (73) 

which gives, with insertion of allowed T and F, all possible constraints on the color combinations in the 
sense of the vector z'm (l44l) . From the non- vanishing possibilities we list only combinations relevant for 
four-quark condensates and contract them to achieve relations between components of z: 

r = l,f = l =^ = -2zi + 2z2 + Zi + 2^6 - 2^8 , 
r = 75, f = 75 ^> = -2zi - 2z2 + Z4- 2zq - 2zs , 

= -2zi + Z2-ZQ + 2z8 , (74) 
= -2zi + 22:2 - 4^3 + Zi- Az5 - 2zq + Azj + 2zs , 

r = 2757", r = 75 =^ o = izg + ziQ. 

This set of constraints is equivalent to (|45l)-(|49l) in section [3TT1 

6.3 Four-Quark Expectation Values in the Nucleon 

Supplementary to Tab. [3] we collect the underlying coefficients to be understood in connection with the 
work of Drukarev et al. [|53l1 . 



r = 2757", r = 2757/3 
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Expectation 
value 


Parameters in [531 

N = p N = n 


Minimal Modification 

N = p N = n 


Mean Value 

N = s^ 


V,uu 

T,uu 

iT.uu 
bN 

A,uu 

jtP,UU 


3.94 4.05 
0.52 0.51 
-0.13 -0.02 
0.98 1.02 
0.05 < 0.01 
-0.45 -0.50 
-0.06 -0.01 
1.91 1.96 


3.939 4.047 
0.520 0.510 
-0.143 -0.023 
0.968 1.009 
0.045 0.007 
-0.471 -0.502 
-0.054 -0.009 
2.002 2.030 


3.993 
0.515 
-0.083 
0.989 
0.026 
-0.487 
-0.032 
2.016 



Table 4: Coefficients of pure flavor nucleon four-quark expectation values (in units of 
{'ll)vac ~ (—0.245 GeV)'^) as determined in |53] in tlie terminology introduced there 
and modified values from a fine-tuned parameter set which fulfill the constraints (l45l) - 
( |49l) . The parameters k^^'^ and k™^'^ are finally derived from the right column which 
shows the result for isospin symmetric baryonic matter 
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Expectation 


Parameters in [53^] 


Mean Value 


value 


N = V 


N = n 


2 


jjSud 


3.19 


3.19 


3.19 


V,ud 

Ql J\T 


-0.44 


-0.44 


-0.44 


iy,ud 

iV 


-0.29 


-0.29 


-0.29 


T,ud 

Qi AT 

iV 


0.19 


0.19 


0.19 


N 


0.18 


0.18 


0.18 


A,ud 

(X AT 

iv 


0.43 


0.43 


0.43 


7 A,ud 


-0.06 


-0.06 


-0.06 


TJ-PyUd 


-0.20 


-0.17 


-0.185 


TjVS,ud 


-0.28 


-0.21 


-0.245 



Table 5: As Tab. |4] but for coefficients of nucleon four-quark expectation values 
parametrizing mixed flavor structures as determined in 115 3n and the mean values used 
to calculate medium strength parameters k^^'^ in isospin symmetric matter The modifi- 
cations referring to pure-flavor four-quark condensates are not needed here. 
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